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Abstract 

We report on a study of the B DD*j{23l7) and B DD*j{2'i60) decays based on 123.7 x 10^ 
BB events collected with the Belle detector at KEKB. The B DD* j{2317) and B DD*j{2A60) 

decays have been observed for the first time. Wc observe the Z)*j(2317) decay to Dgir^ and 
i:>*j(2460) decay to the D*tt^ and Dg-f final states. We also set the 90% CL upper limits for the 
decays L>;_;(2317) ^ L>*7, £»;_;(2460) L>*7, D*j(2460) ^ DsTT^ and L>;_;(2460) DsTT+tt-. 
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I. INTRODUCTION 



Recently the BaBar collaboration reported on the observation of a new Dgn^ resonance 
with a mass of 2317 MeV and a very narrow width A natural interpretation is that this 
is a p-wave cs quark state that is below the DK threshold, which accounts for the small 
width 0. This interpretation is supported by the observation of a D*7i^ resonance Q by 
the CLEO collaboration Both groups observe these states in inclusive e~^e~ processes. 
The mass difference between the two observed states is consistent with the expected hy- 
perfine splitting of the p-wave Ds meson doublet with total light-quark angular momentum 
j = 1/2 01 . However, the masses of these states are considerably below potential model 
expectations and are nearly the same as the corresponding cu states recently measured 
by Belle The low mass values have caused speculation that these states may be more 
exotic than a simple qq meson system 0, S S lit To clarify the nature of these 



states, it is necessary to determine their quantum numbers and decay branching fractions, 
particularly those for radiative decays. In this connection it is useful to search for these 
states, which we refer to as D*j, in exclusive B meson decay processes. 

We search for decays of the type B — > DD*j, which are expected to be the dominant 
exclusive D*j production mechanism in B decays. Because of the known properties of the 
parent B meson, angular analyses of these decays could unambiguously determine the D*j 
quantum numbers. Moreover, since QCD sum rules in HQET predict that p-wave mesons 
with j = 1/2 should be more readily produced in B decays than mesons with j = 3/2 [l3| . 
the observation of B DD*j would provide additional support for the p-wave nature of 
these states as well as serve as a check of these predictions. 

In this Letter we report on a search for the B DD*j{2317) and B DD*j{2A60) 
decays based on a 123.7 x 10^ produced BB pairs at the KEKB asymmetric energy e^e~ 
collider [l^. The inclusion of charge conjugate states is implicit throughout this report. 



II. EVENT SELECTION 

The Belle detector has been described elsewhere l3]- We select well measured charged 
tracks that have impact parameters with respect to interaction point (IP) that are less than 
0.2 cm in the radial direction and less than 2.5 cm along the beam direction (z). We also 
require that the transverse momentum of the tracks be greater than 0.05 GeV/c in order to 
reduce the combinatorial background from low momentum particles. 

For charged particle identification (PID), the combined information from specific ioniza- 
tion in the central drift chamber (dE/dx), time-of-flight scintillation counters and aerogel 
Cerenkov counters is used. Charged kaons are selected with PID criteria that have an effi- 
ciency of 88%, a pion misidentification probability of 8%, and negligible contamination from 
protons. All charged tracks with PID responses consistent with a pion hypothesis that are 
not positively identified as electrons are considered as pion candidates. 

Neutral kaons are reconstructed via the decay Kg — >■ n~^7r~ with no PID requirements for 
the daughter pions. The two-pion invariant mass is required to be within 9 MeV/c^ (~ 3a) 
of the mass and the displacement of the 7c'^tt~ vertex from the IP in the transverse 
(r — (p) plane is required to be between 0.2 cm and 20 cm. The direction in the r — (p plane 
from the IP to the 7T~^7r~ vertex is required to agree within 0.2 radians with the combined 
momentum of the two pions. 

Photon candidates are selected from energy deposit clusters in the Csl electromagnetic 
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calorimeter. Each photon candidate is required to have a laboratory energy greater than 
30 MeV with no associated charged track, and a shower shape that is consistent with an 
electromagnetic shower. A pair of photons with an invariant mass within 12 MeV/ (? (~ 2.5cr) 
of the 7r° mass is considered as a 7r° candidate. 

We reconstruct Dt mesons in the 07r+, and decay channels. The mesons 

are reconstructed from the K^K~ pairs with the invariant mass within 10 MeV/c^ (~ 2.5r) 
from the mass. The K**^ mesons are reconstructed from K~tt^ pairs with an invariant 
mass within 75 MeV/c^ (l-Sr) of the K*^ mass. After calculating the invariant mass of 
the corresponding set of particles, we define a Ds signal region as being within 12 MeV/c^ 
(~ 2.5cr) of the Ds mass. mesons are reconstructed in the D* —>■ Ds'j decay channel. 
The mass difference between D* and Ds candidates is required to be within 8 MeV/c^ of its 
nominal value (~ 2.5a). 

We reconstruct D^{D~) mesons in the i^'+vr", K^-k^-k^-k^ and K^tt^ti^ (K^ti^ti^) decay 
channels and require the invariant mass to be within 12 MeV/c^ of the D^{D^) mass. For 
the 7r° from the 5° K+ 71 7r° decay, we require that the 7r° momentum in the T(45') 
center-of-mass (CM) frame be greater than 0.4 GeV/c in order to reduce combinatorial 
backgrounds. 

We combine D candidates with a D^*^^ and a vr^, 7, or tt+tt^ pair to form B mesons. 
Candidate events are identified by their CM energy difference, AE = Ei) — E^, and the 
beam constrained mass, Mbc = \/E'^ — i^iPiY-, where E^, is the beam energy and Pi and Ei 
are the momenta and energies of the decay products of the B meson in the CM frame. We 
select events with M^c > 5.2 GeV/c^ and \ AE\ < 0.2 GeV, and define a B signal region of 

5.272 GeV/c^ < M^c < 5.288 GeV/c^ and | < 0.03 GeV. In the cases with more than one 

(*)+ 

candidate in an event, the one with the D and Dl masses closest to the nominal values is 
chosen. We use a Monte Carlo (MC) simulation to model the response of the detector and 
determine the efficiency (l6l |. 

III. BACKGROUND SUPPRESSION 

Variables that characterize the event topology are used to suppress background from the 
two-jet-like e"'"e~ qq continuum process. We require | cos6'thr| < 0.80, where 6'thr is the 
angle between the thrust axis of the B candidate and that of the rest of the event; this 
eliminates 77% of the continuum background while retaining 78% of the signal events. To 
suppress combinatorial background we apply a restriction on the invariant mass of the D 
meson and the 77° or 7 from D^j decay: M{Dtt^) > 2.3 GeV/c^ M{Dj) > 2.2 GeV/c^. 

IV. RESULTS OF THE ANALYSIS 

A. Calculation of branching fractions 

The AE and D*j candidate's invariant mass {M{D*j)) distributions for B DD*j 
candidates are presented in Fig. ^ where all D^ and D~ decay modes are combined. Each 
distribution is the projection of the signal region of the other parameter; distributions for 
events in the M{D*j) and AE sidebands are shown as hatched histograms. 

Clear signals are observed for the D D*j {2320)[Ds7i^] and DD*j{2A60)[D*7r°, Ds^f] final 
states. The measured masses for the L';j(2317) and L',*j(2460) are (2319.8±2.1±2.0) MeV/c^ 
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FIG. 1: AE (left) and M{D*j) (right) distributions for the B DD^j candidates: (a) 
L'*j(2320) ^ Ds-K^, (b) Dlj{2m) D*tt° and (c) D*j{2460) Ds'j. Points with errors repre- 
sent the experimental data, hatched histograms show the sidebands and curves are the results of 
the fits. 



and (2459.2 ± 1.6 ± 2.0) MeV/c^. 

For each decay channel, the AE distribution is fitted with a Gaussian signal and a linear 
background function. The Gaussian mean value and width are fixed to the values from a 
MC simulation of signal events. The region AE < —0.07 GeV is excluded from the fit to 
avoid contributions from other B decays of the type B — >■ DD*jX where X denotes an 
additional particle that is not reconstructed. The fit results are given in Table HI where 
the listed efficiencies include intermediate branching fractions. The statistical significance 
of the signal quoted in Table HI is defined as ■>/— 2 ln(£o / J^max) , where Cmax and Cq denote 
the maximum likelihood with the nominal and with zero signal yield, respectively. 

The combined fit quoted in Table |n] is the summed results for the and D~ modes 
assuming isospin invariance. The normalization of the background in each sub-mode is 
allowed to float while the signal yields are required to satisfy the constraint Ni = N^q ■ 
B{B DD*j) ■ Ei , where the branching fraction B{B DD*j) is a fit parameter; A''^^ is 
the number of BB pairs and Ei is the efficiency, which includes all intermediate branching 
fractions. From the two B DD* j{2460) branching fraction measurements, we determine 
the ratio B{D*j{24:60) Ds-f)/B{D*j (2460) D*n°) = 0.38 ± 0.11 ± 0.04. 

The signals for the B DD* j (2317) [DsTI^] and B D D*j {2460) [D* ti^ , Ds-f] channels 
have greater than 5a statistical significance. Figure |21 shows the AE distributions for the 
other channels, where significant signals are not seen. We set 90% confidence level (CL) 
upper limits for these modes based on an event yield N that is calculated from the relation 

C{n)dn = 0.9 C{n)dn, where C{n) is the maximum likelihood with the signal yield 
equal to n. 
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TABLE I: Branching fractions B DD*j decays. 
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TABLE II: Combined fit results. 
Decay channel B, 10^^ Significance 
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B. Angular analysis 

The helicity angle analysis can provide information about spin of the decaying particle. 
We study the helicity angle distribution for the D*j(2460) Dg'y decay. 

The helicity angle Oos^y is defined as the angle between the D*j(2460) momentum in the 
B meson rest frame and the Dg momentum in the D*j(2460) rest frame. The distribution 
shown in Fig.lHlis consistent with MC expectations for a J = 1 hypothesis for the D*j(2460) 
(x^/n.d.f= 5/6), and contradicts the J = 2 hypothesis (x^/n.d.f= 44/6). 

C. Cross checks &; systematic uncertainties 

We study the possible feed-across between all studied D*j decay modes using MC. We 
also analyse a MC sample of generic BB events corresponding to our data sample. No 
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FIG. 2: AE distributions for the decay channels with not significant signals: (a) D-D*j (2320) [D* 7], 
(b) DD*j{2m)[D*'y], (c) LiD*j(2460)[Ds7r+7r-], (d) DD*j{2m)[D,TT°]. Open histograms repre- 
sent the experimental data and curves show the results of the fits. 
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FIG. 3: The D*j(2460) — > Dsj helicity distribution. The points with error bars are the results of 
fits to the AE spectra for experimental events in the corresponding bin. Solid and dashed curves 
are MC predictions for a J = 1 and J = 2 hypothesis, respectively. The highest bin has no events 
because of cut to the invariant mass. 



peaking background is found. 

As a check, we apply a similar procedure to decay chains with the same final states: 
B —>■ DD*, B —>■ D*Ds and B — > D*D*. For each mode, we measure branching fractions 
that are consistent with the world average values [l7| . 

The following sources of systematic errors are considered: tracking efficiency (1-2% per 
track), kaon identification efficiency (1%), vr° efficiency (6%), Kg reconstruction efficiency 
(6%), efficiency for slow pions from D* Dtt decays (8%), D branching fraction uncer- 
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tainties (2%-6%), signal and background shape parameterization (4%) and MC statistics 
(3%). The uncertainty in the tracking efficiency is estimated using partially reconstructed 
D*^ D^iT~^, K^Ti^Ti~ decays. The kaon identification uncertainty is determined 

from D*^ — s> D^Tr"*", —* K^tt~^ decays. The 7r° reconstruction uncertainty is obtained 
using decays to K~^7t~ and K^ti^ti^. We assume equal production rates for B^B~ and 
B^B'^ pairs and do not include the uncertainty related to this assumption in the total sys- 
tematic error. For the calculation of the branching fractions, the errors in the meson 
branching fractions are taken into account. These uncertainties are dominated by the error 
of the Ds (pTT^ branching ratio of 25% The overall systematic uncertainty is 30%. 

V. CONCLUSION 

In summary, we report the first observation of -B DD* j{2320) and B DD*j{2A60) 
decays. The measured branching fractions with the corresponding statistical significances 
are presented in Table ITTl The angular analysis of the D*j(2460) Dg'j decay supports the 
hypothesis that D*j(2460) is a 1+ state. 
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